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. Abstract. The debate over whether kilometer-sized sohds, or planetesi- 

' mals, assemble by collision- induced chemical sticking or by gravity-driven 

unstable modes remains unsettled. In light of recent work showing that 
Q I gravitational growth can occur despite turbulent stirring, we critically 

•7- . evaluate the collisional hypothesis. Impact speeds in protoplanetary disks 

I reach ~ 50 m/s in a laminar disk and may be larger in turbulent disks. 

. We consider the role of elastic and plastic deformations, restructuring 

of "rubble piles," and sticky organic matter. Coagulation is possible for 
^ ■ dust grains with large surface area-to- volume ratios. Simple energetic ar- 

guments show that bouncing, cratering, and fragmentation should domi- 
■ nate collisional dynamics between millimeter and kilometer sizes. 
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^ I 1. Introduction 

Qh' Solid planets and gas giant cores form in a hierarchical growth process spanning 

Q . some 13 orders of magnitude in radius from interstellar dust grains to spectac- 

I ular giant impacts. An important demarcation occurs at kilometer sizes, where 

^ ■ coupling to the gas disk is weak and substantial escape speeds of 1 m/s allow 

growth by inelastic binary collisions. Two theories aim to explain the origin 
of planetesimals from smaller solids with negligible surface gravity and strong 
gas couplings. This paper examines what we will call the Collisional-Chemical 
Sticking Hypothesis (hereafter CCSH), that km-sized bodies arise from a long 
\ chain of agglomerative growth, due to chemical forces during inelastic collisions. 

The gravitational instability (hereafter GI) hypothesis appeals the collective self- 
gravity of a dense midplane layer to collect small solids into larger planetesimals 
(Goldreich &; Ward 1973, Youdin & Shu 2002, hereafter YS). The characteristic 
wavelength of GI predicts planetesimal sizes ^ few km for disks with enough 
mass to produce a solar system like our own. Leapfrogging intermediate sizes by 
direct assembly from ^ cm solids avoids most problems associated with sticking 
and rapid radial drift (see §2.1.). 

The main obstacle to GI is that turbulence opposes the settling required to 
create a dense midplane layer. Even an otherwise laminar disk generates mid- 
plane turbulence, since particle settling introduces destabilizing vertical shear 
(Weidenschilling & Cuzzi 1993, hereafter WC). Sekiya (1998) overcame this ob- 
stacle by appealing to the enhanced buoyancy of stratified, high metallicity disks. 
YS showed that turbulent motions can only suspend a finite surface density of 
solids, 5]p,crit- GI occurs in a dense midplane layer when the actual Sp exceeds 
this saturation threshold. However the required surface density ratio of solids 
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to gas, Sp,crit/5^g; i-e. the disk metallicity, is about an order of magnitude over 
solar abundances. Various enrichment mechanisms were described in YS. A 
particularly robust finding is that aerodynamic radial drift has the global conse- 
quence of enhancing Sp as the particle disk shrinks in radius. Youdin & Chiang 
(2003) explored these "particle pileups" in detail, including inertial effects and 
self-consistently generated turbulent stresses in the particle layer. 

Since GI can overcome the heretofore fundamental obstacle of vertical shear, 
we assess the problems facing the CCSH. YS argue against the plausibility of the 
CCSH with evidence based on micro-gravity experiments, the size distribution 
of meteoritic inclusions, the limited sticking ability of icy planetary ring debris, 
and "intuition." Our goal here is to physically explain the "intuitive" notion 
that collisions between macroscopic solids with weak surface gravity, at speeds 
characteristic of freeway driving, are not cohesive. 

The most fundamental obstacle to sticking is energetic: the entire rel- 
ative kinetic energy of two bodies at large distances, i^oo = fnv^f^, must 
be dissipated, where m is the reduced mass and Voo is the relative speed at 
large distances. This requires some combination of (a) efficient energy loss, 
and (b) binding energies, C/bind < 0, which increase kinetic energy on im- 
pact: Kiiap = + I C/bind I- The dissipative efficiency needed for sticking, 
i^oo/-^imp = 1/(1 + &)■, drops appreciably below unity when the focusing factor 
6 = |C/bind|/-f^oo ^ 1- Gravitational focusing, which scales as ^ oc E? /v1^, affects 
bodies with sufficiently large radius, R. For small bodies, surface energies, such 
as van der Waals interactions, give oc /v^ , where the size falloff ranges from 
w = —1, for a maximal contact area ~ i?^, to w = —5/3, for elastic compression 
(see §3.1.). Over a large size range, extending from mm — 10 m, no relevant 
binding energy exists. Until an incredibly efficient, yet non-destructive, dissipa- 
tion mechanism is identified, the CCSH will lack a solid physical foundation. 

This proceeding is organized as follows: §2. investigates collisional speeds 
expected from laminar flow (§2.1.) and turbulent fluctuations (§2.2.). Section 
3. explores collisional binding mechanisms including elastic (§3.1.) and plastic 
(§3.2.) deformation as well as restructuring within "rubble piles" and the role of 
organic matter (§3.3.). Section 4. contains concluding remarks. 

2. Collision Speeds 

Before planetesimals form, gravitational scattering is negligible. Particle random 

velocities arise from frictional coupling to the gas, including kicks from turbulent 
fluctuations, drift induced by the gas disk's sub-Keplerian rotation, and vertical 
settling, which becomes small near the midplane as vertical gravity vanishes. 

2.1. Laminar Drift 

Aerodynamic gas drag gives particles radial, v^-, and azimuthal, v^, deviations 
from Keplerian circular speeds, vk = (Nakagawa, Sekiya, & Hayashi 1986): 



Obstacles to Collisional Planetesimal Growth 



3 




Pg/P 

Figure 1. Contours of in-plane relative velocity induced by laminar 
gas drag, labeled in units of r]VKp/Pg ~ 54/3/pg m/s, between particles 
with stopping times Tstop.i and rstop,2- 



where r is the disk radius, J7 the Kcplcrian frequency, r] = —dP/dr/{2pgQ^r) 
measures centrifugal pressure (P) support; p = pp + pg \s the spatial mass 
density of solid particles, pp, plus gas, pg] and Tstop = ^^^stop measures the 
coupling of solids to gas, with Tgtop 1 1) implying strong (weak) coupling. 
The stopping time, istop; depends on the drag law (e.g. Epstein or Stokes) that 
applies for given particle sizes and gas properties (Weidenschilling 1977). 

The velocity scale, rjVK ~ 54 m/s, is independent of r in the minimum 
solar nebula model (see YS), hereafter MSN. Different temperature profiles, 
T oc r~'^, give rjVK oc r^l'^~i. Our discussion ignores the inertial scaling factor 
Pgl p < 1, which is 0(1) unless pp > pg, in which case GI should develop (YS). 
Particles with Tstop ~ 1 experience the fastest radial drift, Vr,max ~ —V'^K- 
These marginally coupled particles have a size Rri ~ 1 m, at 1 AU in the MSN, 
giving the well-known and problematic drift timescale: ~ AU/(r/'Ux) ~ 100 yr. 
Elsewhere a meter is less significant. At larger r, pg decreases giving weaker 
coupling and smaller Rti values. Beyond about 3 AU, Tstop = 1 particles lie 
in the Epstein regime (size less than the gas mean free path) and obey Rt\ = 
Sg/(\/27rps) ~ 10 cm(r/10 AU)^'^/^. Note that as any particle tries to grow 
past the local Rri value, to slow its radial drift, it continually moves to regions 
where Rri is even larger. This struggle could be more futile, if shorter lived, 
than Sisyphus's since all motion is downhill (radially inward). 

Figure 1 (generalized from figure 3 of WC) plots the relative velocity, arising 
from Vr and drift, between particles with a range of stopping times, and thus 
sizes. Only similarly sized particles avoid 0{r]VK) relative motions. Impacts 
above a threshold speed, Vcr will not give coagulation. While Vcr uncertain and 
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size dependent, Vcr <S 50 m/s is reasonable (as discussed in §3.). Thus colli- 
sional growth requires maintaining a narrow size distribution. This is a difficult 
proposition since many processes; e.g. coUisional cascades, radial migration, and 
vertical settling, broaden the size dispersion. 

2.2. Turbulent Fluctuations 

Protostellar disks must be at least intermittently turbulent for stars to accrete. 
This turbulence is often modeled as an anomalous viscosity, = ac^/U, where 
Cg is the gaseous sound speed and the uncertain viscosity parameter takes values 
10~^ < a < 10~^ to reproduce T Tauri accretion rates. Solids couple to available 
turbulent eddies, augmenting the particulate velocity dispersion, Cp. 

Assume a Kolmogorov spectrum whose largest eddies have a turnover time 
to ~ and thus a typical speed vq ^JTm^ ^ 30y^a3(r/AU)~"^/^ m/s, where 
otz = a/lO"^. Particles with tstop/^o ^ '^stop < 1 acquire Cp ~ fo relative to the 
gas. Larger particles, with Tstop 3> 1, couple poorly to all eddies, and proceed 
with Cp <C vq on nearly Keplerian orbits. The dispersion of particles relative to 
each other is as large as vq only when one of the particles is marginally coupled, 
Tstop ~ 1 (or when the two particles have Tstop,! <C 1 <C Tstop,2)- Pairs of small 
particles, both with Tgtop 1, remain so tightly coupled to eddies that they 
have little motion relative to each other (Markiewicz, Mizuno, &: Volk 1991). 

Turbulence, unlike laminar drift, imparts relative motion to particles with 
identical stopping times. These motions are the largest, ~ vq, near Tstop ~ 1- 
This closes the narrow "bridge" of low speed collisions in Figure 1 (also see 
Figure 4 in WC), which we already noted would be difficult to cross by coUisional 
growth. The difficulties facing the CCSH escalate further in a turbulent disk. 

Small grains, on the other hand, can coagulate in turbulent disk, especially 
those with tgtop less than the eddy turnover time at the dissipation scale, to = 
to/y'Reo, where Rcq is Reynolds number for the largest eddies. This regime, 
where the flow is laminar and Brownian motion determines Cp, corresponds to 
sizes, 

R < (27r)-V4^y^ « 100(r/AU)-3/4a3 /xm, (3) 
where a is the atomic cross section and /x is the mean molecular weight. 

3. Contact Mechanics 

To evaluate the CCSH, we need to understand the aftermath of collisions at 
speeds up to ~ 50 m/s. Do they result in coagulation (mass gain), fragmen- 
tation (mass loss), or restitution (no change in mass)? We consider identical 
spheres with radius R and density ps, unless otherwise noted. CoUisional out- 
comes depend on material properties, the most relevant are probably Young's 
modulus, E; the yield strength, Y; and the binding energy per unit area, 
7. Other influential factors include impact parameter, non-sphericity, surface 
roughness, and viscoelastic creep, and are considered only briefly. 

At low impact speeds, materials deform clastically, i.e. rcvcrsibly. In this 
regime, the ratio of applied stress (force per unit area) to induced strain (relative 
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compression) remains constant at E. When imposed stresses exceed Y, deforma- 
tion becomes plastic, or permanent, resulting for instance in dents. Elastic and 
plastic deformation during impacts generates a contact area, A ^ ttR'^. Chemi- 
cal forces, mainly van der Waals interactions, generate a surface binding energy, 
t^bind = —7^) where 7 varies from 30 — 300 — 3000 ergs/cm^ for quartz, ice, 
and iron grains, respectively. 

We argued in §1. that sticking is problematic (requires very efficient energy 
dissipation) for 9 < I. Assuming a maximal contact area, A ~ ttR^, the upper 
limit on 9 due to surface forces is: 



Koo Ps,3 R Vm/s 

with ps,?, = Ps/(3 g cm~^) and 7300 = 7/(300 erg cm~^). This suggests that 
coagulation between cm-sized solids at speeds above 10 cm/s is difficult. More- 
over, surface roughness reduces effective binding forces by factors ^ 3 because 
only crests of asperities come in contact (Johnson 1985). 

3.1. Elastic Collisions 

Elastic sticking is an oxymoron if one defines an elastic collision as having a 
coefficient of restitution, e = 1. However, kinetic energy is lost during elastic 
deformation since elastic waves (EWs) absorb kinetic energy (-ftTimp) during com- 
pression, but do not perfectly return deformation energy during rebound. These 
wayward EWs, which one can think of as wandering from the impact region, 
eventually dissipate, though not necessarily during the abrupt collision. 

An elastic theory of coagulation was developed by Chokshi, Tielens & Hol- 
lenbach 1993 (hereafter CTH) who were mainly interested in the coagulation of 
small (^ mm) grains for which surface binding energies are significant. They 
calculate a critical approach speed, Wcr, above which sticking will not occur: 

^5/6 f R \ ^^/^ cm 



E^/^R^/^pl''^ Vmm 

Equation (5) is physically plausible since v^r increases with surface binding en- 
ergy, 7; decreases with R, as the surface area-to- volume ratio falls; and decreases 
with Young's modulus, E, since stiffer materials share smaller contact areas and 
place less energy into elastic waves. The threshold speed is quite low. Equation 
(5) gives numerical values for ice, while speeds for graphite and quartz are low- 
ered by factors of 6.5 and 30, respectively. Thus elastic effects cannot explain 
growth beyond mm-sizes which have drift speeds v,. ~ l(r/AU)^/^ cm/s. 

Nevertheless, an order-of-magnitude "derivation" of (5) is instructive. Press- 
ing a flattened area. A, on a surface with radius of curvature, R, requires a 
displacement, 5 ~ A/R, perpendicular to the surface. Young's modulus relates 
stresses, F/A, imposed by a static force, F to the local strain, ~ 5/\/~A, giv- 
ing A ~ {FR/Ef/^ and 8 ~ [F"^ / {RE'^)^/^ . We treat compression statically 
because elastic wave periods are short compared to the coUisional duration. 

The force that interests us is the surface binding force, which is equal and 
opposite to the force needed to separate two spheres, Fggp = —Sir'-fR (Johnson, 
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1985, p. 127). This result, surprisingly independent of E, follows from the 
fact that ?7bind induces a finite contact area between spheres in the absence of 
external forces. Equating this surface energy, —jA, to the work required to 
separate the bodies, 

We can also estimate equilibrium values for 6 ~ {■y'^R/E'^y/^ , A ^ {'fR? / E)"^/^ , 
and, most importantly, C/bind ~ -'-^^/^R'^I^E''^/^ . 

Our heuristic argument that sticking needs K^o <, | C/bind I (unless dissi- 
pation is very efficient, which should not be the case for clastic deformation) 
reproduces equation (5) to order unity. CTH achieved their result differently. 
They estimated the amount of energy lost to EWs, C/iost, by numerically evolv- 
ing a spectrum of EWs in response to time-varying collisional stresses. They 
found that U\ost is a fixed fraction of the binding energy, ?7iost ~ -^iC/bindl! and 
Koo < U\ost gives equation (5). 

The best hope for the CCSH probably involves large bodies sweeping up 
small grains, which dominate the relative kinetic energy and binding area. How- 
ever, the large relative velocity between unequal bodies, r]VK ~ 54 m/s, implies 
via equation (5) that only extremely small "grains" could stick: 30 — 3 — 0.5 nm 
for ice, graphite, and quartz grains, respectively. However, impact speeds might 
be reduced below 54 m/s as the gas flow is deflected around the larger body. 
Nevertheless, maintaining a reservoir of small sticky grains requires destructive 
collisions which limit or even reverse the growth of the "runaway" body. So it 
remains to be shown whether this growth mechanism has the proper sign. 

3.2. Plastic Deformation 

Car crashes show that plastic deformation can dissipate large amounts of kinetic 
energy in violent collisions. This may not be directly relevant to planetesimal 
formation. Only sufficiently ductile materials can permanently change shape in 
response to large stresses, others simply rupture. The former category includes 
metallic elements and alloys, while the latter includes glass, bone, and crucially, 
rock. Of course planetesimals must form prior to segregation of an iron-rich 
core, and evenly dispersed metallic elements probably cannot sustain plastic 
flow. Would you feel safe driving a car that was 1/4 metal and 3/4 rock? 

Even if planet forming material does deform plastically, we should be aware 
that the regime of stability is rather narrow. Contact pressures only slightly 
above Y give plastic flow that is confined by similarly large elastic stresses to a 
small region. Somewhat larger stresses, ~ 2>Y ^ give "uncontained" plastic flow 
that breaks through the surface allowing displaced material to escape (Johnson 
1985). During planetesimal formation, uncontained flow can cause mass loss by 
sub-sonic cratering and, at high enough energies, shattering, though more study 
of these outcomes is deserved. Since confined plastic fiow requires a narrow 
regime of contact pressures, its relevance to the CCSH, which exposes solids to 
orders of magnitude variations in impact stresses, is not clear. 

Plastic fiow develops at an impact velocity, vy, calculable from elastic Hertz 
theory (which ignores surface binding) by equating i^oo to the elastic deforma- 
tion energy, / Fd5. Requiring the maximum stress to be Y gives : 
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or ~ 10 cm/s for medium soft steel. Thus most collisions between metals do 
indeed involve plastic deformation (Johnson 1985, p. 361). 

3.3. Aggregate Restructuring and Viscoelastic Organics 

So far, we have implicitly assumed that collisions involve uniform, compact ma- 
terials. What if sticking relies on the detailed structure of impactors? For 
instance porosity (i.e. "fluffiness") should enhance dissipation. However open 
structures are fragile, and moreover would be compacted by repeated impacts. 
More plausibly, growing solids might be loosely bound "rubble piles," i.e. an ag- 
gregate composed of many smaller units. Restructuring of the sub-units during 
impacts could be quite dissipative. By analogy, collisions between bean bags 
(sacks filled with beans or foam pellets) have a lower coefficient of restitution 
than collisions between individual large beans or pellets. Of course bean bags 
have a cloth cover that astrophysical rubble piles lack. 

Blum (1990) constructs a model along these lines by assigning low values of 
E = 2 X 10^ dyn/cm^ and Y = 100 dyn/cm^ to readily deformable aggregates. 
A serious problem exists with models of this type. Impact velocities, VcoU > 
vy ~ 10^'^ cm/s exceed the yield stress, and the sound speed in this material 
is only yj Ej « 2 m/s. Impacts in the cm/s — 10 m/s range would be highly 
disruptive events, inducing "uncontained" plastic flow and possibly driving shock 
waves into the material. Blum ignores these possibilities and extrapolates results 
(e.g. for e) from the moderate speed regime, p^v^^^jY <^ 10~^, all the way to the 

hypervelocity regime, Ps^cou/^ ^ i where they do not hold (Johnson 1985, 
p. 366). 

These caveats aside, Blum's model gives a maximum sticking speed 0.5 < 
Vcc < 5 m/s for his mm-sized aggregates, depending on impact parameter. The 
predicted coagulation velocity limit for ices is quite low, Vct ~ 3(i?/mm)"^ mm/s 
(as read from the graph). Amazingly (since the included physics is quite differ- 
ent) this agrees very with equation (5) as evaluated for ice, including the drop 
of Vci with R. Blum developed this model to explain experimental results which 
showed that mm-sized wax coated spheres (with a glass/ wax mass ratio of 2 — 
5) have a 20 it 12% sticking probability for Voo ~ 1 m/s. It is a mystery how 
such results (low sticking probabilities, small sizes and moderate speeds) com- 
bined with the theoretically expected decrease in Va with R could be interpreted 
as supportive of the CCSH. The sizable wax fraction invites skepticism, but it 
justified by the existence of organics in meteorite falls. 

Along these lines, Kouchi et al. (2002) collided 5 mm copper spheres into 
targets with a mm-thick coating of organics, including PAHs and oxygen- and 
nitrogen-containing materials. They report temperature-dependent sticking at 
speeds up to Vcr ~ 5 m/s at 250 K. The threshold fell to ^ 2 m/s and 
^ 1 m/s at 300 K and 200 K, respectively. Kouchi et al. offer viscoelastic 
creep, a time-dependent relation between stress and strain, as an explanation 
for enhanced stickiness. The role of plastic deformation in the copper is unclear. 
The observed coagulation speeds are probably insufficient to make planetesimals, 
in particular due to the expected decrease of Vcr with size. Furthermore, the 
experimental conditions are not likely to be repeated in protostellar disks. The 
targets had mm-thick layers of pure organics, while carbonaceous chondrites 
(the most carbon rich meteorites) have a bulk carbon content of only ~ 3.5%. 
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4. Conclusion 

We studied the implications of particle velocity dispersions in protoplanetary 
disks for the collisional growth of planctesimals. Over a wide range of sizes, mm 
— 10 m, sticking is strongly inhibited. Available binding energies are negligible 
compared to kinetic energies, and a suitably efficient, yet non-destructive, dis- 
sipation mechanism has not been identified. The most elegant solution to this 
problem is the Goldreich-Ward hypothesis, which uses collective gravity to avoid 
this problematic coagulation regime. 

Particle concentration mechanisms, such as vortex trapping, secular insta- 
bilities, and quiescent nodes in turbulent flow (see YS for details) may play a 
role in planetesimal formation. Higher densities favor gravitational instability, 
and could also promote collisional growth if collisional speeds arc significantly 
reduced. Regarding the possibility of long-lived particle-trapping vortices, only 
marginally coupled particles (rgtop ~ 1) are efficiently trapped. Smaller particles 
remain attached to streamlines and larger bodies can plow through the vortex. 
This limits the overall reduction in relative velocities. 

For many years, supporters of the Goldreich-Ward hypothesis were on the 
defensive due to valid concerns about the stirring provided by turbulence. Pro- 
posed solutions to this problem (Sekiya 1998, YS) should reestablish the viability 
of particulate gravitational instabilities. To make the debate more productive, 
coagulation proponents should similarly acknowledge, and explain how to over- 
come, the considerable energetic obstacles to sticking. 
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